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Abstract 

We achieve a detailed understanding of the n-sided planar Poisson- 
Voronoi cell in the limit of large n. Let Pn be the probability for a 
cell to have n sides. We construct the asymptotic expansion of log pn 
up to terms that vanish as n — > oo. We obtain the statistics of the 
lengths of the perimeter segments and of the angles between adjoining 
segments: to leading order as n — > cxd, and after appropriate scaling, 
these become independent random variables whose laws we determine; 
and to next order in 1/n they have nontrivial long range correlations 
whose expressions we provide. The n-sided cell tends towards a circle 
of radius (n/47rA)2, where A is the cell density; hence Lewis' law for 
the average area An of the n-sided cell behaves as An ^ cn/\ with 
c = \. For n — >• oo the cell perimeter, expressed as a function R{(j)) of 
the polar angle (j), satisfies (PR/dcj)'^ = F{(j)), where F is known Gaus- 
sian noise; we deduce from it the probability law for the perimeter's 
long wavelength deviations from circularity. Many other quantities 
related to the asymptotic cell shape become accessible to calculation. 
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1 Introduction 



Experimental data on cellular structures are often compared to the properties 
of a Voronoi diagram jl], sometimes also called a Voronoi tessellation. This 
well-defined mathematical object is a partitioning of space into convex cells 
constructed around 'point particles' (also called 'centers,' 'seeds,' or 'germs') 
in such a way that each point of space is in the cell of the particle to which it 
is closest. When the distribution of the point particles is uniformly random, 
the result is said to be a Poisson-Voronoi diagram; it is to naturally occurring 
tessellations what the ideal gas is to real gases. 

In this work we will be interested in planar cellular structures. An ex- 
ample analyzed more than seventy years ago by Lewis ^ is the epidermal 
epithelium of cucumber, which is composed of many flat, tile-like cells. For 
other applications of Voronoi diagrams outside of physics we refer to the 
encyclopedic work by Okabe et al. [3]. 

Within physics, two-dimensional cellular structures appear in many dif- 
ferent experimental situations, of which we mention only a few. Cerisier et al. 
studied the cellular patterns in surface-tension driven Benard convection. 
Elias et al. |S] studied two-dimensional magnetic liquid froth. Sire and Seul 
jHj analyzed droplet patterns in late-stage coarsening of a two-dimensional 
binary mixture. Moriarty et al. [T analyzed nanostructured cellular lay- 
ers by means of the Voronoi construction. Earnshaw and Robinson [HI E] 
performed a Voronoi cell analysis of the structure of two-dimensional col- 
loidal aggregates. Much important experimental work on tessellations comes 
from two-dimensional soap froths, as studied experimentally by Glazier et al. 
PH] and by Stavans and Glazier In simulations of the two-dimensional 
solid-liquid transition (recent examples are Zahn et al. ^2] and Quinn and 
Goree [12]), Voronoi cells constructed around the particle positions define 
the nearest-neighbor relations between atoms, which in turn allow the iden- 
tification of lattice defects. Richard et al. 14^ used the Voronoi construction 
for the geometrical characterization of hard sphere systems. 

Evidently, there are many different reasons why none of the above systems 
should be described exactly by a Poisson-Voronoi diagram. In dense fluids the 
particle positions are correlated and the resulting Voronoi diagrams are not of 
the Poisson type. Both soap froths and binary mixtures coarsen with time, 
and the experimentally determined properties of these cellular structures 
usually refer to a stationary scale invariant distribution (sometimes called 
the 'equilibrium distribution') observed at long times. Clearly such systems 
are best analyzed, in principle, by a dynamical theory such as developed, e.g., 
by Marder ^Sl, by Stavans et al. 16,, and by Flyvbjerg |17]. Nevertheless, 
comparison of dynamical experimental data to the static Poisson-Voronoi 
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diagrams may give a clue as to what elements should enter into a correct 
dynamical theory. 

Discrepancies between the experimental data and the predictions of a 
Poisson-Voronoi diagram are often interpreted as the consequence of some 
physical, biological (see e.g. ^3), or other action that one would like to 
identify. It is therefore of importance that the properties of Poisson-Voronoi 
diagrams be understood as well as possible. 

In field theory 'pure' Voronoi diagrams play a role since attempts were 
made pHlEOl to construct a theory on a lattice of randomly located sites. The 
Voronoi construction then again determines the nearest neighbor pairs among 
the sites. Such a random lattice has the advantage of being statistically 
invariant under arbitrary translations and rotations while still preserving a 
short- distance cutoff. Field theory provides still another reason for interest 
in Voronoi diagrams: Godreche et al. [21] showed that the ensemble of cells 
with trivalent vertices in a plane is related to the problem of counting planar 
Feynman diagrams ('Feynman foam') of a field theory with cubic interaction. 

1.1 Statistics of planar Poisson-Voronoi cells 

Planar Poisson-Voronoi diagrams are characterized by their statistical prop- 
erties. The first and foremost of these, whether in experimental studies, in 
theoretical work , or in simulations, is the probability p„ that a cell have n 
sides. Others include the average number {n) of cell sides; the average area 
An of an n-sided cell; the statistics of the cell perimeter and of its angles; and 
correlations between neighboring or nearby cells. The topology of the plane 
imposes that (n) = 6. One may readily obtain several of the other statistical 
properties analytically by expressing them as integrals on the positions of 
the point particles [221 El- The calculation of the fraction p„, however, has 
still defied solution [201 E]- The distribution p„ peaks at n = 6 and falls of 
rapidly for large n. For n = 3 Hayen and Quine have expressed as 
a fivefold integral which they evaluated numerically to seven decimal places. 
For n = 4,5,... the values for pn obtained from Monte Carlo simulations 
have been tabulated. All Monte Carlo studies cited in reference jB] (see e.g. 
[23 123 I2n]) were restricted to n 14, which is the largest number of sides 
observed in simulations by conventional algorithms: one has pi4 ~ 10~^. 

As early as 1984 Drouffe and Itzykson [27], by a then novel method, 
estimated the probability pn for n up to 50 by Monte Carlo simulations. 
This work is still today the main reference for the large n behavior. The 
presently known values of pn have their uncertainty typically in the third 
digit for 4 < n < 10 and in the second digit for 10 < n < 16. Very large 
cells with n values higher than 60 were observed by Lauritsen et al. [2H1 
in simulations of a 'Voronoi cell Hamiltonian' which, although of uncertain 
physical relevance, is interesting because it favors the appearance many-sided 
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cells; this procedure does not, however, provide estimates for p„. 

In view of its importance as a system of reference, we study in this work 
the planar Poisson-Voronoi diagram. We obtain a great many new analytic 
results for n-sided cells in an asymptotic expansion for large n. 

Expansions around idealized models are common in physics and as such 
our approach needs no specific justification. The fact that n = oo is not 
physically accessible in any realization known of today, is of no concern, 
since the expansion goes into the physical domain n < oo. After the work is 
done, we will point out the relevance of our results for finite n. In particular, 
the vast body of experimental and numerical results has led to two empirical 
laws that are satisfied to a good extent by many systems: Lewis' law j2] and 
Aboav's law ^j. The first one will be discussed at the end of this paper; the 
second one will be investigated separately |30j as an extension of our present 
work. 

An announcement of part of our results, including a concise indication of 
their derivation, has appeared earlier Here we present the full calculation 
together with many more results. This is Part I of a series of two papers. It 
presents our analytic work and the results that derive immediately from it. 
In Part II we will develop arguments which, with the benefit of hindsight, 
will help to understand our results heuristically and to employ them in other 
applications. 

1.2 Method 

We consider a square L x L subdomain of the plane M^. In this domain 
point particles are placed at random positions, chosen independently and 
with uniform probability. We consider the limit A^, L — oo with the particle 
density N/L"^ = A fixed. In mathematical terminology this is the Poisson 
point process of intensity A in the plane. The density A may be trivially 
scaled away. The Voronoi cells of this point particle system are necessarily 
convex polygons; their number of sides n may be any integer > 3. 

We ask about the statistical properties of the Voronoi cell of an arbitrarily 
selected particle, whose position is taken as the origin of the coordinate 
system. Let the other particles, indexed by a label a, have position vectors 
2Ra. The perpendicular bisectors of these vectors pass through the midpoints 
Ra, and the line segments that constitute the cell sides lie on these bisectors. 
Out of the full set of bisectors only a finite number v that pass sufficiently 
close to the origin will actually contribute a segment to the perimeter of 
the cell containing the origin. We denote these 'contributing' midpoints as 
Ri, R2, . . . , Rjy. They determine the positions of the cell vertices, which 
we indicate by Si, S2, . . . , Sj,. The angles between two successive midpoint 
vectors Rm-i and R^ [vertex vectors S^-i and S^] will be denoted as 
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[as rjm]- The probability Pn that we have i/ = n for a given n, i.e., that 
the selected Voronoi cell have exactly n sides, can be expressed as a 2n- 
dimensional integral on Ri, R2, . . . , R„. 

We evaluate this integral for asymptotically large n by a multidimensional 
steepest descent method. Although the principle of this approach is simple, 
its execution is beset by numerous interdependent complications. The neces- 
sity to deal with these is the reason for the length of this paper. We were 
initially guided by the idea of describing the perimeter in terms of a Markov 
process as a function of the polar angle; the Markov aspect of our description 
will, however, stay in the background and we will not seek to elaborate on 
it. 

We write p„ = Tre""^, where EI is a 'Hamiltonian' of purely geometrical 
origin and Tr stands for the integrations. For a suitably chosen 'free' or 
'noninteracting' Hamiltonian Hq we set H = Hq + V and write 

p„=(Tre-»o)(e-^)„, (1.1) 

where (. . .)g denotes the average with respect to e"'^". We will refer to the 

two factors, Tre~''^° and {e~^)^, as the 'noninteracting' (or 'free') and the 
'interacting' problem, respectively. For clarity of presentation, rather than 
for technical reasons, wc study them separately. 

In the course of our analysis the appropriate Hq appears to be a sum of 
noninteracting terms in the angles and rjm' for n — > 00 these variables 
decouple. The evaluation of Tre""^" amounts to a 2n-dimensional steepest 
descent problem in the positive orthant {i.e. the subspace with all C,m,Vm > 
0), with contributions concentrated in a boundary layer of width of order 
n~^. The answer can be obtained exactly up to corrections that vanish 
exponentially as n — > 00. 

The result obtained for Tre"''^" is meaningful only if we are able to show 
that the calculation of (e"^)^ actually leads to a controlled perturbation series 
in negative powers of n. The expression for the 'interaction' V is so opaque 
that a priori it is not clear whether it lends itself to any kind of large n 
expansion, let alone how it will contribute to the n dependence of the final 
result. 

We overcome this problem by considering the perimeter segments as a set 
of vectors S^, — S^-i which, by analogy to the theory of elasticity, we expand 
about their values in a regular n-sided polygon. An analysis of considerable 
difficulty shows that we may write V = Vi + n~2 V2, where Vi is a quadratic 
form in terms of the 'elastic deviations', and n~2 V2 collects together a great 
number of higher order terms. The maximum of e~^'^ is located in the 
positive orthant, inside the boundary layer identified in the free problem. 
The 'elastic' interaction Vi has a coupling constant proportional to n^^, but 
is long-ranged along the perimeter in the same way as a one-dimensional 
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Coulomb interaction. There is a full separation of spatial scales: Hq is a 
Hamiltonian of independent 'microscopic' variables; and Vi has its weight 
concentrated on an effectively finite number of wavelengths at the scale of 
the cell. 

The integrations in the interacting problem involve the confluence of the 
boundary layer maximum and the maximum of e""^^ . They can be dealt with 
via a Hubbard-Stratonovich transformation. Due to the long-rangedness of 
the interaction and the weakness of the coupling, an expansion in negative 
powers of n becomes possible. In the end we find that the leading order 
term in the interacting problem is equal to a finite constant: {^~'^\ ^ C as 
n — >• cx). 

Various technical obstacles surge in the course of our analysis. To con- 
clude this subsection, we mention three of them. 

First, in previous approaches [23 IHH IHlj the domain of integration in 
phase space was always defined in terms of a set of n inequalities. Here 
we choose variables of integration such that it becomes fully explicit. After 
suitable transformation the integration bears on only one radial variable and 
2n — 1 angles. The radial integration can be carried out exactly and the 
remaining angular integrations constitute the true problem. There is no 
single set of angular variables in terms of which the integrand takes a simple 
form. 

Secondly, the angular variables satisfy sum rules stemming from various 
different origins. It is indispensable to strictly keep track of each of them. 
To that end we denote them as 'geometrical sum rules,' a 'gauge condition,' 
and a 'no-spiral' constraint; their meanings will be explained as they are 
encountered. 

Thirdly, the asymptotic expansion for n — »• oo requires, as usually, that 
we postulate at the outset the appropriate scaling with n for the variables 
of integration. But since n is also the typical number of terms in many 
summations that occur, our scaling analysis has to deal with these sums, too. 
Some remarkable cases of canceling leading orders need particular attention. 

Throughout this work we use only those methods of mathematical physics 
that ordinarily lead to exact results. 

1.3 Results 

We obtain a detailed understanding of the statistical properties of the n-sided 
Voronoi cell when n is large. The main qualitative results are summarized 
as follows. 

(i) The probabilty p„ for a Voronoi cell to have n sides behaves asymp- 
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totically as 



P„ = -,^[l + 0(n-)], n^oo, (1.2) 



47r2 {2n) 
where 

1 4 
- + - 



c=n(i 

q=l 



= 0.344 347... (1.3) 

and where 0{n'^) denotes a term which for n — > oo is to leading order pro- 
portional to n'^. Alternatively this may be written as 

logp„ = -2nlogn + nlog(27rV) - | log {2^n^C-^n) +0{n'^). (1.4) 

In fll.2j) the factors (47r^)^^(87r^)"/ {2n)\ and C originate from the free and the 
interacting problem, respectively, i.e., from the factors Tre~^" and (e"'^)^ in 
fll.ll) . The constant C in carries a great deal of meaning. Its factor with 
index q stems from the 'elastic' modes with wavelength 27r/g. Hence the 
contributions to C are seen to come essentially from the large scale, long 
wavelength modes. With the decrease of scale they damp out as ~ 1/q'^. 

The remaining results, (ii)-(vii) below, arise as byproducts of the calcu- 
lation of pn. They have, nevertheless, an independent interest equal to that 
of Pn, if not greater. 

(ii) As shown in figure^ we write C,m for the angle by which the perimeter 
turns at the mth vertex (it is also the angle between two successive midpoint 
vectors); and rim for the angle by which it advances as it passes from the 
mth to the (m + l)th vertex. The are positive and add up to 27r, and 
the r]m similarly. Because of these sum rules they have a weak 'background' 
correlation which in the limit n ^ oo vanishes as ~ (here and henceforth, 
the sign ~ denotes asymptotic proportionality). In that limit the and rj^ 
become independent random variables with the identically distributed 
according to the law 

uiO = ^ exp (- ^) (1.5) 



and the rjm according to 



= ^exp (- — ). (1.6) 



These distributions have the averages ^ = f] = 2?™^^, as had to be the case. 
The background correlations have contributions not only from the sum rules, 
but also from the interaction V, which renders them nontrivial. We find the 
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Figure 1: Heavy line: the perimeter of the Voronoi ceh of the particle in the 
origin. Dashed and dotted lines connect the origin to the midpoints and vertices, 
respectively. Right angles have been marked. The figure defines the angles and 
r]m discussed in this section as well as Pm and j^n to be introduced in section 2. 

expressions for these correlations in sections [7| and |H1 but will not state the 
results here. 

(iii) Let the mean value of the midpoint distances i?i,i?2, ■■■,Rn be 
denoted by i?av = J2m=i 

-Rav = Rc + SRav, Rc = (^;^) • (^■'^) 

In the limit n — *• cxo we have (-Rav) = Rc and 5-Rav has the Gaussian distri- 
bution 

p((5i?av) = (8A)^exp(-87rA<50. (1.8) 

i.e., its width is independent of n. 

(iv) The fluctuations of the Rm around their average i?av have the n 
independent variance 

((i?^-i?av)')^-^, n^oo, (1.9) 

with Ca^ = - + 4g-^)-^ = 0.333 418.... Here and henceforth 

(...) indicates an average over all Voronoi diagrams and ~ stands for asymp- 
totic equality. Equations (jl.7p and (jl.9|) together tell us that the Rm all scale 
as n2 and that the line segments linking the midpoints run mostly inside 
an imaginary tube of width ~ along a circle of radius Rc = (47rA)~2n2. 
The cell perimeter is the polygon having not the R^, but the as vertices; 
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however, we show that the distance between the two polygons vanishes as 
~ n~2. Therefore the cell perimeter approaches a circle as n — oo. 

This approach to circularity when the number of sides tends to infinity 
is distinct from other approaches to circularity, such as occur for cells whose 
area tends to infinity (see Hug, Reitzner, and Schneider ^j), or whose in- 
radius {i.e., the radius of the largest inscribed circle centered at the origin) 
tends to infinity (see Calka [32]; Calka and Schreiber |SZ1)- In Part II j32] 
we compare the two types of limits in greater detail. 

(v) Let 27rm?2~^ = and Rm = R{4>)- In the limit n — oo the variable (j) 
becomes continuous and may be identified, up to corrections that vanish as 
n~2, with the polar angle of the perimeter. We set 

R{(j)) = Rc + 6R{(f)). (1.10) 

and refer to 6R{(f)) as the 'excess midpoint distance.' On scales large com- 
pared to the angular scale 27m^^ of the individual perimeter segments the 
following holds. The effect of the Poisson point process on the perimeter 
R{(f)) may be represented by a stochastic second order differential equation 
for (5i?(0), 

^ = F(0), < < 27r, (1.11) 
in which -F(0) is zero average Gaussian noise of autocorrelation 

(F(0)F(0')) = 12A k0 - 00 - ^-^(f-^') ] , (1.12) 

where r(0) is a 27r-periodic function whose integral on [0, 27r] vanishes. The 
solution of (|1.11|) is subject to a periodicity condition and an integral condi- 
tion, viz. 

5R{0) =SR{27r), ^[ d0(5i?(0) = (5i?av, (1.13) 

Jo 

respectively, with 5_Rav distributed according to ()1.8|) . Hence in the n — oo 
limit SR{(f)) is governed by Gaussian statistics. The correlation function 
p.l2|) determines the probability law of the noise F{(f)) and hence, by equa- 
tions ()1.1H) . p.l3p . and ()1.8p . of 6R{(j)). The exphcit result is given in section 



(vi) According to Lewis' empirical law |2l Ej the average area An of an 
n-sided cell grows linearly with n. For Poisson- Voronoi cells, in spite of 
attempts, this linearity with n has not so far been proved [SHllSni- Now the 
results (iii) and (iv) above imply directly that 

An ~ vri?^ = oo. (1.14) 

Hence we have derived here from first principles that Lewis' law holds in the 
asymptotic regime and has coefficient \. 
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2 Writing pn as a 2n-fold angular integral 

In this section we start from the basic 2n-fold integral for p„ and subject it 
to a sequence of transformations that will cast it in a form amenable to an 
asymptotic large n expansion. 

2.1 Basic integral for 

An n-sided cell is completely defined by n midpoint positions Ri, R2, . . . , Rn 
(see figure IH). These are uniformly distributed with density 4A, and Pn may 
be written directly as a 2n-fold integral on their coordinates [211 1201 E3 Ell , 



(4A) 

Pn ■ 



fJdR,... dR„ • • • , Rn) e-^^-^(^--'^"), (2.1) 



n 

where the functions x ^ind A to be defined shortly enforce two conditions 
that must be satisfied if the n points are to define a valid Voronoi cell. In 
order to state these conditions we use the polar coordinate representation 
Rm = {Rm, ^m)- In the integrand of (|2.1|) we may set one of the angles, say 
to zero if we compensate by an extra factor 2tt, and we may order the 
remaining angles according to 

= < $1 < $2 < . . . < $n-l < 271 (2.2) 

if a compensating factor (n — 1)! is introduced. We will denote the n angular 
differences between successive midpoint vectors as 

Cm = $m-$m-l, m = 1, . . . , n - 1, 

e„ = 27r - <l>„^i . (2.3) 

The geometry constrains the to the interval (0, vr) and they must satisfy 
the geometrical sum rule 

n 

J2U = 2n. (2.4) 

m=l 

We find it convenient henceforth to scale the distances such that 

4A = 1 (2.5) 

during the calculations, and restore this factor in the final results. Sym- 
metrizing in the angular difference variables by introducing a Dirac delta 
function we obtain from (j2.H) 

Pn=— / RldR^...RndRn / d^i . . . d^ 5(^1 + • • • + " Stt) X 6^-^. (2.6) 



n 



^0 
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We will also need the positions Sm = {Sm, ^m) of the cell vertices (see figure 
These are auxiliary variables in the sense that they can be expressed in 
terms of the C,e and Re. We are now ready to state the two conditions. 

Condition 1. The first condition is that for a 7^ l,...,n none of the 
bisectors passing through the midpoints Rq have an intersection with the 
Voronoi cell. This condition is fulfilled if these midpoints stay outside of a 
region of the plane equal to the union of n disks centered at and of 
radius \Sm [40j. If this region has area A, then the probability for it to be 
empty of midpoints (other than the n that determine the cell perimeter) is 
equal to exp(— ^) times its area. An explicit expression for A in terms of the 
Rm and may be found by straightforward geometry (see Calka [331 El|)- 
Condition 2. The second condition defines the domain of integration. It 
says that each of the n midpoints Rm should have its corresponding per- 
pendicular bisector contribute a nonzero segment to the perimeter. This 
condition was shown by Calka jSSI IMI to read explicitly 

i?m-isin^m+i + -Rm+isin^m > -Rmsin(^m + ^m+i), m = l,...,n, (2.7) 

with the conventions Rq = Rn, Ri = Rn+i, and ,^1 = ^n+i- The function x 
in fl2.1|) and ()2.6|) is the indicator of the domain in phase space where ()2.7p 
is satisfied. 

This fully defines the integral ()2.6|) for p„. In being a sum on n equivalent 
interacting objects, viz. the midpoint vectors R^, this integral is analogous 
to a partition function in statistical mechanics. But although interesting, 
this observation provides no direct key as to how to go about its evaluation. 

2.2 Angular variables 

There is no single set of variables in which the calculations of this work take 
a simple form. We will, throughout, employ simultaneously different sets of 
angular variables, to be introduced in this subsection. 

The polar angles "^rn of the vertex vectors satisfy an ordering analo- 
gous to (Q, 

< ^1 < ^2 < • • • < < 27r. (2.8) 
There is, however, no particular ordering relation between the two sets 
{"^m \ m = 1, . . . ,n} and {$m | m = 1, . . . , n}. We will not use the {"^m I n^' = 
1, . . . , n} any more, but instead, by analogy to equation (j2.3p . introduce the 
angular differences 

Vm = m = 2,...,n, 

7]i = ^i-^„ + 27r, (2.9) 

which are constrained to (0, vr) and satisfy the geometrical sum rule 

n 

J2vm = 2n. (2.10) 

m=l 
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Figure 2: A configuration of midpoints and Rm+i leading to a negative angle 
7m- The relation /3m+7m = Cm remains valid. The negativity of 7^ is accompanied 
by a counterclockwise outward spiraling of the perimeter. The legend is the same 
as for figure 1. 

Next, as shown in figure [H the set {Pm, 1m\m = 1, . . . ,n} will represent the 
angles between a vertex vector and a midpoint vector. They are defined in 
terms of the and by 

$m„i , m = 2,...,n, 

^m, m = l,...,n. (2.11) 

The Pm and 7^ are geometrically constrained to the interval (— |, |), and 
hence, contrary to the C,m and 7]^, niay be negative. An example of a negative 
7m is shown in figure El Whereas this potential negativity may at this point 
not seem worthy of much attention, it will be seen in section El to have 
important consequences for the scaling of these angles with n. Because of 
the geometrical sum rule 

n 

5^(/5m + 7m) =27r, (2.12) 

m=l 

the set {Pm, ImlfTT' = 1, . . . ,n} contains one redundant angle. Hence, leaving 
out 7„, the set of 2n — 1 variables 

{/3i,7i, . . . ,/5„_i,7„_i,/3„} (2.13) 



Pm = ^m- 

A = ^1, 

7™. 
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is complete in the sense that, given $o = 0, all the other angles $m and \E'm 
are easily expressed in terms of the angles of the set ()2.13|) . 

The set {^m, Vm\fn = 1, • • • , n}, due to the sum rules ()2.4|) and ()2.10|) . 
contains two redundant variables and is not complete: it specifies only the 
angles of the between themselves and of the between themselves, but 
not the relative angle between the two vector systems. To fix their relative 
orientation, one other angle is needed, for which we will take Pi. A second 
complete set of 2n — 1 variables is therefore 

{ei,T/i,...,en-i,r/„-i,/5i}. (2.14) 

The sets (I2.13j) and ()2.14|1 have /3i in common. The transformation relation 
between their 2n — 2 remaining variables is 

^m=Pm + 1m, Vm = + Pm+l , 171 = 1, . . . , U - 1, (2.15) 

and its nonlocal inverse 

m—1 

e=i 

m—1 

7m = -Pi + '^{^e - Vi) + ^m, m = 1, . . . ,n - 1. (2.16) 



where by definition Yl'i=i — summand. It is easy to verify that this 

transformation has unit Jacobian. Relations ()2.15|1 are valid also for m = n 
with the convention that Pn+i = Pi'-, and relations ()2.1(jj) are valid for m = 1 
and m = n, respectively. 

2.3 Making the domain of integration explicit 
2.3.1 Rewriting the indicator function 

From the elementary observation that 

Rm-l — SmCOS Pm Rm = >S'm COS 7m (2-17) 

we deduce the relation 

Rm cos 7„ 



im 



Rm-l COS/?„ 



m = l,...,n. (2.18) 



Hence, when Ri, . . . ,Rn are given, Pm determines ■jm and vice versa. We now 
turn to the inequality ()2.7|) and divide it by Rm sin ^ 

mSiTi^rn+i- The ratios 

of the midpoint distances that then appear may be eliminated: we express 
them in terms of the angular variables by means of relation (|2.18p . After 
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some trigonometry and ^m+i cancel out of the inequality, which takes the 
form tan 7m + tan/3m+i > for m = 1, 2, . . . , n. Since — | < 7m,/?m+i < 
|, it follows that the much simpler and geometrically necessary positivity 
condition 

7m + /5m+i>0, m = l,2,...,n, (2.19) 

is equivalent to condition (j2.7j) . Therefore (j2.19j) defines the domain of inte- 
gration in phase space. 



2.3.2 Transforming from to (3m and 7^ 

Next, considering the outer integration variables . . . , in ()2.6|) as simple 
parameters, we will transform the integrals on the angles in two succes- 
sive steps. First, we consider 7^ as a function of j3m defined by (|2.18|) . so 
that in view of ()2.15b ) we have = Pm + arccos[(i?m/-Rm-i) cos /5m]. We 
transform from to the new integration variable jSm, taking due account of 
the Jacobian = d^m/^^Pm , which is equal to 

U = "'•''f"''^"' . (2.20) 

We order the (3^ integrations such that, going from the outermost one inward, 
we encounter successively those on We integrate Pi on its full 

interval (— |, |). For given Pi equation (j2.18|) fixes 71. Inside the Pi integral, 
in order to satisfy condition (j2.19p for m = 1, we integrate P2 from —71 to |. 
For given P2 equation ()2.18|) fixes 72, and so on. In this way the conditions 
fl2.19|) are all satisfied, except the one at the end which requires 7^ + /3i > 
and which we will therefore incorporate below by means of a Heaviside theta 
function. After all this equation (j2.6|) becomes 

Pn = — RidRi . . . RndRn / d^ / d/?2 • • • / d/3„ 

'0 J-7T/2 J-y. -'-7„_i 



n 



9i^n + Pl)jlj2...Jn5{J2iPm + lm)-2n) C""^, (2.21) 



X ^ 

m=l 



in which 7^,, wherever it occurs, is a function of Pm given by ()2.18p . and 
where the indicator x is no longer needed. The passage from ()2.6|) to ()2.2ip 
has broken the symmetry between the Pm and the 7^ and introduced an 
orientation in the cell perimeter. We have moreover broken the cyclic invari- 
ance in the index m by imposing a fixed order of integration of the angular 
variables. These symmetry breakings first appeared in an initial approach 
IH] where we constructed the perimeter segments by a Markov process with 
transition probabilities between successive vertices. We will not elaborate on 
this here. 
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We wish to introduce in ()2.21|) the 71, 72, . . . , 7n as additional independent 
variables of integration. One easily verifies that 

sin 7^ / cos 7^ \ 

/ '^7m —Rm-lO{Rm —Rm-l]=l- (2.22) 

COS fJm V COSfJrn ^ 

We now take the product on m of the LHS of this expression and insert it in 
the integrand of ()2.21|) . Still defining 



_ sm{p^ + 7^) 

COS^ (3rri 



we find that ()2.21|1 becomes 



o_ /'OO 

Pn = — RidRi . . . i?ndi?„ / d/5d7 



f n TmRm-iS{R^ - ^^Rm-i)] e^^ (2.24) 

L -'- -'- ^ COS iJm -I 



COS /5, 

m=l 

with our earlier convention Rq = Rn and where we employ the shorthand 
notation 

dA / d7i / d/32 / d72... / d/3„ / d7n 

-7r/2 J-7j J~I3^ J-lr.-! -^-Pn 

n 

X ^(7n + A) 5( 5^(/5,„ + 7m) - 27r). (2.25) 

m=l 

The operator / dpd'j represents the integral on all possible shapes of the 
Voronoi cell without any reference to its radial scale. The fact that the do- 
main of integration appears here fully explicitly constitutes the main achieve- 
ment of this subsection. 

An expression for A in terms of the variables of integration Rm, Pm, and 
7m may be obtained directly from its geometrical definition as the area of 
the union of n disks. After the necessary rewriting one finds 

n 

^=zY.^'rn Him) + a{Pm+l)] , (2.26) 
m=l 

where a is defined by 

a{uj) = uj + tantu + tutan^ u, (2.27) 

which is equivalent, and close in form, to the expression used by Drouffe and 
Itzykson ^27j in their simulations. 
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2.4 Radial integrations 

2.4.1 Introducing the mean midpoint distance i?av 

The average value of the n midpoint distances i?2, ■ ■ ■ ,Rn will be denoted 

by 

n 

i?av = Yl ^rn ■ (2.28) 

m=l 

The circle of radius i?av will act as a circle of reference and we proceed on 
the hypothesis that the area difference A — irRl^ will be small in a sense still 
to be determined |4^. Let us now normalize the Rm according to 



Rm — RavPrn ; (2.29) 

SO that because of ()2.28j) the 'reduced midpoint distances' pm satisfy 



n 

— 1 ^ 

n 

m=l 



'J2Pm = ^- (2.30) 



Instead of ()2.29|) other choices of normalization would have been possible. 
Therefore we will henceforth refer to ()2.30p and equivalent relations as the 
'gauge' condition; this in order to distinguish it from other sum rules that 
play a role. 

We insert a factor 5(-Rav — J2m=i in the integrand of the RHS of 
fl2.24|) and apply an extra integration on i?av to it. The radial integrations 
then take the form 



poo roc POO _ "' _ 

dRi...dRn = / d/^av / dRi... dRn 5{R^^ - n"^ ^ R, 
Jo Jo Jo 

POO POO 

= / di?avi?L"' / dpi...dp„5(l-n-i 

(2.31) 

and result ()2.24|) becomes 
27r 



n 



Pn = — 

n 



POO poo 
/ di^avi^av"' / dpi...dpn5{l-n-^Y.P^) 

Jo Jo 

/n n 
mi \X\KPm~ '^Pm-l)] \ n pIT^ (2-32) 
L COS Pyyi -I L -I 



where by convention po = Pn and in which the angular integrations are 
defined by ^TI^ . 
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2.4.2 Integrating on the Rm at fixed i?av 

The n — 1 integrations on pi, . . . , p^-i may be carried out successively and, 
due to the delta functions in the integrand, cause p^, wherever it occurs, to 
be replaced by 

cos7mCOS7m_i . . .COS71 

Pm = n n ^ Pn , 771=1,...,^-!. (2.33) 

COS Pm COS Pm-l ■ ■ ■ COS Pi 

There remains one delta function inside the product on m in ()2.32|) . We 
rewrite it as 

/ COS71 . . .cos7„ \ 5{G) 

•^i^" R irP^) = o — ' (2-^^) 

V cos pi . . . cos fjn / 27rp„ 



where 



1 

G = — ^ ( log cos 7„ - log cos (3m) . (2.35) 



27r 

m=l 



The integration on p„ in ()2.32|) . finally, is easily performed due to the delta 
function constraint in the first line of that equation and leads to the cancel- 
lation of the factor in ()2.34|1 . The result is that ()2.32j] becomes 

Pn = - / d/3d7 5(G) [ n p^tJ rdR^^Rl:-' e-^, (2.36) 

in which now the reduced midpoint distances pi,...,p„ should be viewed 
as functions of the integration variables determined by the n — 1 equations 
(IZTIHll together with condition ((2301) • 

It appears from ()2.36|) that the 2n angles of integration satisfy a nontrivial 
supplementary constraint G = 0. This constraint is easily interpreted: for 
an arbitrary set of angles (/?,7) = {/3m, 7m I = the perimeter 

will spiral instead of close onto itself after a full turn of 2tt; the factor S{G) 
enforces its closure. We will refer to it as the 'no-spiral' constraint. We see 
that equation ()2.34|) imposes that ()2.33|) hold for m = n. Hence requiring 
the validity of ()2.33|) for all m = 1, . . . , is another way of incorporating the 
no- spiral constraint. 

In what follows we will extensively use an auxiliary variable represent- 
ing the excess length of the reduced midpoint distances pm over their mean 
value, 

Pm = l + Tm, m = l,...,n. (2.37) 
From (j2.3(J|) one deduces that they satisfy 

n 

$^r^ = 0, (2.38) 

m=l 

which is an alternative expression for the gauge condition. 
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2.4.3 Integrating on the mean midpoint distance i?av 

Equations ()2.29p . ()2.30p . and ()2.33|) allow us to express Rm as i?av times an 
m dependent function of the angles which, if desired, may be made explicit. 
It is then clear from (|2.26p that A is Rl^ times a function only of the angles. 
Bearing in mind that we expect A to be close to the area vri?^^ enclosed by 
the circle of reference, and anticipating later developments, we write 

A = TTRl^{l+n-W), (2.39) 

which defines V. Integrating on i?av then converts (j2.36j) into 

Pn = ^^^^ /"d/3d7 5(G) [Hp^tJ [77(1+ n-V)]-" (2.40) 

m=l 



The integrations that remain in the RHS of ()2.40|) bear exclusively on the 
angular variables, i.e. on the shape of a cell with average midpoint distance 
scaled to unity. 

For completeness and for use in section El we now determine V explicitly. 
Substituting (jTTTTfl in and subsequently (IT^ and (pTTfl in (jT^ . 

comparing to ()2.39|) . and still using the geometrical sum rule f|2.12|) . we obtain 



V{(3,'y) = — ^(l + 7:„)^[tan7^-7„ + tan^„+i 



m=l 



n 

n 



+ -fm tan^ 7m + Prn+l taU^ Pra+l] 

+ ^ + ^-)(^- + (2.41) 

m=l 

We note that ()2.4H) is antisymmetric under a change of sign of all the angles. 
2.4.4 The probability distribution of i?av 

In this subsection we sidestep and make an observation that will be ex- 
ploited only in section 1^1 Although the outcome (|2.4U|) of the integral on i?av 
is known, something more may be learned by applying a steepest descent 
analysis to it. Assuming, as is suggested by the notation n~^V in ()2.39p and 
as will be confirmed in section 15.3.21 that V{l3,'-f) remains of order when 
n ^ oo, we then find that the integral on i?av draws its main contribution 
from a Gaussian centered at the saddle point -Rav.c and having a variance 
f^av,c given by 

<c(A7) = [4vr(l + n-V)]-i = i^(l + 0(n-i)). (2.43) 
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The results fj2.42|) and ()2.43p are to leading order in n independent of (/?, 7). 
The 0{n~^) corrections in these equations do depend on (/3, 7) but will play 
no role in this work. Hence, provided we confirm the scaling of V^, we have 
found here that the n-sided Voronoi cell has a characteristic radial scale of 
order na Furthermore, if we define 

= T^-\vh^ R^^ = R^{1 + n-3rav), (2.44) 

then it follows from the above that in the limit n — 00 the variable r^v has 
the probability distribution 

Pav(rav) = exp {-2rl) + 0{n--^). (2.45) 

With 5i?av = n~^Rcrs_^ and with the midpoint density 4A restored this be- 
comes equation (jl.8|) of the introduction. 



2.5 Passing to the angles of integration and 77^ 

So far the f3m and 7^, [the set ()2.13p ] have been the variables most convenient 
to work with. However, in order to prepare for the large n expansion, it is now 
necessary that we pass to the set of variables (|2.14|) consisting of the and 
rjm as well as of Employing the notation (^, rj) = {^rn, Vm\fn = ^,---,n}, 
^ = {^m \ ^ = 1, • • • , n}, and r] = {rjm \m = 1, . . . , n}, we now proceed as 
follows. 

By writing G{^,r];f3i) we indicate explicitly that we want to consider G 
as a function of ^, rj, and f3i. The prime on G' stands for differentiation with 
respect to Pi. Let furthermore P*{^,ri) denote the solution of the no-spiral 
constraint 

G(e,r7;/3.)=0. (2.46) 

In ()2.40p one easily converts the integration on (/5, 7), defined in ()2.25|) . to one 
on the new variables (^, rj) by shifting the integration intervals appropriately 
(except the outermost one, for jSi, which remains unchanged) and employing 
the fact that 7^ -|- /5i = 27r — X]m=i Vm- With the aid of the identity 

e{2n-J2vm)= dr]n6{J2vm-2n) (2.47) 

m=l m=l 

and using the occasion to rearrange some factors in the integrand we recast 
()2.40|) in the form 



in 

Pn = — 



X 



/ dA d^dr^ 5( V - 27r) 5( V r^^ - 27r) 
||^[nplT,„]wl^„-Vr" (2.48) 
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where 



T/2+/3, A-7r/2+7-^ rn/2+t3^ 

d^dr] = I d^i / dT^i / d6 . . . 



^0 ^0 

dr/„_i / de„ / d7]n. (2.49) 

^0 JO 

The notation is hybrid; the variables 71, /?2, 73, • • • , /3n in the upper integration 

hmits in ()2.49|) should be viewed as functions of ^, 77, and /3i given by ()2.16|) . 

Each limit of integration depends only on the variables of more outward 

integrations, as of course it should be. 

We note that it is not possible at this stage to carry out the integral on 
since /?* depends on ^ and rj and the integration limits of the integrals 

/ d^drj in turn depend on [3i. 

2.6 Free and interacting problem 

2.6.1 Defining a free problem and an interaction 

In order to relax the constraints ()2.4|) and ()2.10j) enforced by the first two 
delta functions in ()2.48p . we will exploit the integral representation 5{X) = 
J^^^^ds exp(— sX), valid for any constant c. Upon introducing a parameter 
K which at this point is arbitrary and scaling conveniently with n we can 
write 

= ^(/tSm^m + (1 ^ f^)^mVm ~ 27r) 5 (Tim^m ~ ^mVm) (2.50) 
2 pc+ioo pioo 



-27ri/ 

where we choose c > 0. Substituting ()2.50|) in ()2.48|) we find 

2 pc+ioo pioo /"""/S 

' / ds dt 

' c—ioo J ~ioQ J—tt/2 



/ ri rc+ioo pica pTT/A n 

Pn = [TT^) ds dt d/?i e2-"Mdedr/5(/?i-/?.)e-=(^'*) 

V27ri/ Jc-ioo J-ioo J-7r/2 J 

= Tre-'^^^'*), (2.51) 



in which the last line defines the Tr operator, and the 'Hamiltonian' HI is 
given as the sum of a 'free' Hamiltonian Ho and an 'interaction' V, 

m{s,t) = Mo{s,t) + V, (2.52) 
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where 



g-Ho(s,t) _ ^-n 1^ ■J~|' g-ns[<m+(l-K)r;m]~nt(Cm~»?m) (2 53) 
m=l 

n 

e-^ = \G\^,r^-p^)\-'[\[plT^i;^]{l + n-^Vr\ (2.54) 

m=l 

The factor inserted in p.53|l and compensated in p.54|l serves to ensure 
that T^^^ has a finite hmit as 0. The Hamiltonian Hq defined by 

fl2.53|l is obviously is sum of 2n independent single-variable terms. However, 
the reason why the splitting ()2.52p should be exactly as defined by ()2.53p and 
fl2.54|) . and not something else, will have to be borne out by the calculation. 

With equations ()2.51|) - ()2.54|) we have achieved the purpose of this section. 
They have nothing of the simple elegance of other model Hamiltonians in 
statistical mechanics and are at first sight hardly suggestive of any particular 
analytical approach. However, the work of this section has served to prepare 
the ground for what follows. 

So far our calculation is valid for all finite n = 3,4,.... The initial 
problem ()2.H1 is exactly equivalent to problem ()2.5H1 . However, we cannot 
continue beyond the present point: it is impossible to carry out the integral 
on Pi in ()2.51|) since depends on the and rjm whose upper integration 
limits in turn depend on f3i. This entanglement will be lifted in section El as 
soon as we start considering the limit n — oo. 

2.6.2 Splitting the calculation 

Ensuant on equation (j2.52|) we have the formal factorization 

Tre-==(Tre-=o)(e-^)o, (2.55) 
whence follows the splitting 

logp„ = logpW + \og{e-\, (2.56) 

where we employ the alternative notation pn^ = Tr e~'^°. We will refer to the 
calculation of pn^ as the 'free' problem and to that of (e"'^)^ as the 'interacting 
problem.' Although technically not strictly necessary, it will greatly enhance 
the clarity of presentation to treat these problems in two successive steps. 

In the course of our work we will also be interested in correlation functions 
of the angles. We write (F)^ for the 'noninteracting' average of any function 
F of the angular coordinates, i.e., 

TrFe-=» 
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A full average (F) may be reexpressed in terms of noninteracting averages 
according to 

In section 0] we will solve the noninteracting problem in the limit n ^ oo. 
The study of the interaction V will be taken up in section [Ul It will require 
considerable work to show that the second term in ()2.56p . log(e~'^)g, has a 
well-controlled expansion in powers of n~^. Correlation functions will be 
considered in sections [7| and |H1 



3 Scaling 

3.1 Scaling of and r]m for large n 

In the preceding section, where n was an arbitrary integer, we have estab- 
lished the starting point for an expansion for asymptotically large n. We will 
now make the first steps in this expansion. A major problem is to know in 
advance how the various quantities will scale with n. For — > cxd we expect 
all angles to become small and we begin by introducing 

^m = n~^Xm, r]m = n~^ym, m = l,...,n. (3.1) 

This scaling of and rjm is certainly expected, since they constitute two 
sets of n positive variables adding up to 2n. 

3.2 Scaling the integral for j9„ 

The integral fd^dr] included in the trace in ()2.51|) is defined in ()2.49|) . It 
there appears that the upper integration limits of the C,m and rjm, given in 
terms of the Pe and 7^, are all positive. Hence after passing to the scaled vari- 
ables Xm and Hm we may replace these upper integration limits by 00. Since 
the integrals converge exponentially, the error should be expected to vanish 
exponentially with n. It then becomes possible to commute the integration 
on Pi through those on the Xm and ym and carry it out with the aid of the 
factor S{Pi — P^). Henceforth the angle Pi, wherever it still occurs (such as in 
expressions for V to be given later), should be viewed as the function P^{C,, rf) 
given by the solution of (j2.46j) . 
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As a result the integral ()2.51|) takes the simplified form 



c+ioo 



SKXj-fL -5(1 ■ l^^y-m' t{^m ym} 



m=l 771=1 







X 



exp(-V(x,i/)), (3.2) 



where we indicated explicitly that V, given by ()2.54|1 . depends on the Xm and 
i/m- We use the ~ sign to indicate that this is an asymptotic equality. 

Provided that V be sufficiently well-behaved, the integrals on the Xm and 
Hm will converge at the condition that 

Re(Ks + t)>0, Re{{l- k)s -t) > 0. (3.3) 

Since Ret = and Res = c > 0, this condition is satisfied if we choose k, in 
the interval < k < 1. It is now of great advantage to transform to variables 
of integration u and v defined by 

U = KS + t, v={l — K)s—t. (3.4) 

In terms of these equation (j3.2p becomes 

V27ri/ 2n 



X 



/"OO POO 

m=l m=l 



X exp(-V(x,y)), (3.5) 

where u and v run independently from — ioo to ioo while passing to the 
right of the origin. The difficulty of evaluating (j3.5p is now hidden in the 
expression for V. 

The final remarks of this section concern practical procedure. Since ()3.5p 
differs from ()2.51|) only by terms that are exponentially small in n and since 
we are preparing for a perturbation series in inverse powers of n, the scaled 
integral (j3.5p together with expression (j2.54|) for V is as good a starting point 
as the original integral ()2.51|) together with ()2.52|) - ()2.54p . In particular, (i) 
setting V = in in ()3.5p yields pn^; and (ii) the average (F)^ defined in 
()2.57|) is, up to exponentially small differences, equal to expression ()3.5|) with 
F inserted and V set to zero, and divided by pn^ . 
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4 Free problem 



4.1 Calculation of 

We will solve the noninteracting problem by calculating p'n' starting from 
expression ()3.5p in which we set V = 0. The 2n integrals on the Xm and ym, 
as well as those on u and f , then factorize. After carrying out the former we 
find 

2n \2-Ki J ; V27ri J 



„ . n^" 27r 27r "-^ 

2n (2n-l)! (n - 1)! 

2\n 



57r 



47r2(2n)!' 



(4.1) 



where in the first step we did the integrations by closing their contours around 
the poles in the origin. Additive corrections to (|4.1|) should be expected to 
vanish exponentially with n. 

It is instructive and useful for later to analyze the integrals also by the 
steepest descent method. In the case of the u integral the integrand may be 
written as exp(— 2nlogu + 27rnM), which has a saddle point for u = Uc = tt"^. 
Upon taking its second derivative, it appears that contributions come from a 
neighborhood of extension ~ n~2 around the saddle point. A similar analysis 
holds for the v integral, for which the saddle point is at f = f c = (27r)~^. 
The steepest descent analysis allows us to conclude that any function h{u, v) 
independent of n and inserted in the integrand contributes an extra factor 
h{uc, Vc) to the value of the integral. This rule will be applied in section ED 



4.2 Averages and correlations 

The averages of the scaled variables Xm and i/m are necessarily in any order 
equal to 

(xm) =x = 2n, {yra) = y = 2%. (4.2) 

One obtains the noninteracting results for arbitrary multivariate moments of 
the Xm and ym by starting from ()3.5|) with the appropriate insertions, setting 
again V = 0, and going through calculational steps similar to those above. 
We state the results only for the second moments, 

127r^n 2 Svr^n 

2n + 1 ^ n + 1 



^^-)o = TT—T, ivDo = z—^- (4-3) 



Due to the symmetry of the noninteracting expression [equation ()3.5|) with 
V = 0], the two-point correlations (xiXm) and {yiym) must be index inde- 
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pendent when i ^ m. This observation combined with the geometrical sum 
rules (Q and (jmUl yields 

{xtx^\ = {yeym)o = iT^m. (4.4) 

Furthermore {x£ym\ = {xi^iym)^ for all I and m. Let us define 

5Xm = Xm-X, 6y,rn = Vm - V- (4.5) 

From what precedes we then have, still expanding (j4.3j) in powers of n^^, 

{6xe6xm), = 27c^{6,^-n-'){l + 0{n-')), 

{SyiSyml = 47r2(5,^-n-i)(l + 0(n-i)), 

{Sx,5yJ^ = 0. (4.6) 

These variances obey the geometrical sum rules ()2.4|) and ()2.1U|) . But al- 
though they include terms of order n~^, we will see in sections [7| and |H1 that 
the interaction V generates additional order contributions. 



4.3 Probability distribution of and 77^ 

The probability for the angle to be between ^ and ^ + d^ is given by 
{^{^m ~ We set ^ = n~^x and Fourier represent according to 

1 r°° 

5{xm-x) = — d/i e^^(^^-^). (4.7) 
27r y_oo 

We insert this in ()H.5|1 . commute the /i integral to the outside, and pass 
through the same steps that led from ()H.5j) to the first line of 1)4.11) . The 
result is similar to ()4.ip but with an insertion v?{u — i/i)~^, 

- x))„ ^ - / d/ie-'^- (4.8) 



2n V27ri y u^"" ^ u - i^i' ) \2-ni J / 



Because of our discussion on insertions at the end of section 14.11 the second 
line in ()4.8|1 is equal to pn\uc/ {uc — i/u))^- Upon dividing by pn^ and setting 
Mc = 1/vr we obtain the probability distribution of x in the limit n — > 00, 

{6{Xrr, - X)) = — d/i 



271 J_^ (1 - i7r/i)2 
4exp(--). (4.9) 
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In terms of the original variable ^ = n^^x this becomes the distribution u{^) 
announced in equation (jl.Sp of the introduction. In a fully analogous way an 
insertion 5{rjrn — rf) produces 

(%m-y))o = ^exp(-|-), (4.10) 
from which v{ri) of equation (jl.fij) follows. 

5 Further scaling 

5.1 Scaling of (3^ and 7^ for large n 

The noninteracting calculations of section |3] required no scaling other than 
the one postulated in (jHUl), which concerns the ^.^ and rjm- We will show 
now how the scaling of all further variables follows from this initial postulate 
plus what we have learned about the independence of the and rjm- 

The first one of equations ()2.16p expresses (3m — /3i as a sum of 2m — 2 
terms ^£ and 7]^ which, in the noninteracting case, are all independent. Since 
m is typically of order n, it follows that [3^ — Pi scales as n~2 . But since Pi 
is in no way statistically different from the other Prn, they must all scale as 
n~2. In a similar way one deduces from the second one of equations (|2.16|) 
that the 7^ scale with n~2. We therefore set 

Pm = n~hm, 'ym = n~^Cm, m = l,...,n, (5.1) 

as well as 

P*i^,v) = n~h^{x,y). (5.2) 

Hence the scaling factor for the angles C,m and rim is different from the one 
for the angles Pm and 7^. By symmetry and because of the geometrical sum 
rule ()2.12p we must have (Pm) = (7m) = vr/n. The difference between the 
scaling of these averages and of the typical values of Pm and 7^, which are 
~ ^"2, can arise only because the Pm and 7^ may take negative values, as 
illustrated in figure |21 

It will be useful to have at hand the scaled equivalent of ()2.16|) . Rendering 
explicit the fact that here Pi = /9*(^, rj) we find 

m— 1 

bm = b^{x,y) -n'^^{xe-yi), (5.3) 

m— 1 

Cm = -b^{x,y) +n'^2^(xe - ye) +n'^Xm, (5.4) 

e=i 

valid for m = 1, . . . ,n. 
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The two sums YlT=i -^^ ^"^^ Xl^i Vt- be viewed as the displacements 
of two independent random walks, each built up out of independent incre- 
ments. Figure ^ shows that these displacements describe the progression of 
the polar angles of the midpoint vector and the vertex vector S^, re- 
spectively. Hence the cell perimeter may be viewed as the result of a 
construction based on two entangled random walks. 

5.2 Cancellations of orders in 

The nontriviality of this scaling picture comes from the fact that various fre- 
quently occurring combinations of the angular variables exhibit cancellations. 
For later use we collect some of them here, employing a hybrid notation in 
which on the one hand the 6^ and c^, and on the other hand the Xm and Um 
appear. First of all there are the sums of odd powers of the type 

7m + /3m+l = Vmn'^, 

ll + = ilm + f3m+l) [Ull + Pl+l) - |(7m + 

= \ym,{.C^ + ^m+l)^ ^ ~ iz/m^ i (5-5) 

and analogously 

f^l + ll = \^m{hl + cl;)n-^-\xln-\ (5.6) 
Secondly, we will encounter differences of even powers, 

= Xi{bi-Ci)n~^, 

= -2xebiCe{be - c^)n~^ + - ce)n~i . (5.7) 

These relations are easily extended, if needed, to sums and differences of 
higher powers of the angles. 

5.3 Solution of r/; = 

The solution /3*(^, r/) of the no-spiral constraint G{^, rj; /?*) = has not been 
needed explicitly. Since f5^ plays an important role, we will show in this 
subsection how it may be obtained in the scaling limit. The derivation also 
establishes its uniqueness within the large n perturbation expansion. Finally, 
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because of equivalence by symmetry, the study of (3^ = (3i gives us access to 
the statistics of all angles (3rn and 7m. 
From ()2.35p we have straightforwardly 

m=l 

n 

= ^J2K-cl) +Oin^l), (5.8) 

m=l 

The leading order term is ~ n~^. Its scaling, as well as that of the correction 
term, may be determined with the aid of equations (|5.7|) and taking into 
account the sum on m. For convenience of notation we rewrite ()5.4|) just for 
this occasion as hm = K + bm and Cm = — + Cm, substitute it in the RHS of 
fl5.8|) . set G = 0, and solve for b^:{x, y) using the geometrical sum rule ()2.12|) . 
The result is that 

1 " 

^ n m—1 n 

4™ 2 m=l 1=1 m=l 

where in the second line we returned to standard notation. Inside the square 
brackets of (|5.9p the double sum is of order na and the single sum of order n. 
By taking the average (. . on both sides of (j5.9|) and using (|4.2j) - (j4.4j) one 

finds, after cancellation of the leading order terms, that (6*)^ = vm"^. This 
relation should in fact hold to all orders in n and leads to {Pm)^ = (7m)o = 
nn"^, as had to be the case. The distribution of around its average is on a 
scale which is larger by a factor 122. Restricting ourselves now to this larger 
scale, we substitute Xm = x + 5xm in ()5.9|) and, reasoning as before, keep 
only the leading order. The result, 

n 

b^x, y) = -n-2 ^ m{6x^ - 6ym) + 0(n"^), (5.10) 

m=l 

shows that b^:{x,y) is a weighted sum of n weakly dependent random vari- 
ables and therefore, in the limit of large n, has a Gaussian distribution. 
Its noninteracting variance may be calculated from ()5.10p and ()4.6|) and is 
given by (6^)^ ~ ^vr^. However, the interaction V yields an additional (9(ra°) 
contribution to this variance. Its explicit value to order may be found 
analytically with the aid of the interacting theory of the next sections, but 
we will not pursue this calculation. 
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6 Interaction 



The total Hamiltonian H was split in equation ()2.52|) into a noninteracting 
term Hq and an interaction V, with an ensuing splitting for logp„ in ()2.5(i|l . 
In section 0] we dealt with the noninteracting problem and obtained, in equa- 
tion ()4.H) . the first term of ()2.56|) . This result is useful only if we can also 
control the second term in ()2.56|) . which requires the evaluation of {e'^)^. 
Furthermore, if we wish to study correlations, {F e'^)^ is needed. We will 
now consider these quantities. 

6.1 General 

The interaction V is defined by expression ()2.54|) . which is particularly opaque. 
We now assert that V it is in fact of order We will show this by con- 
structing an analytically tractable expression Vi such that (4^ 

V = Vi+n~^V2 (6.1) 

with Vi of order n° and n~~^Y2 representing a collection of higher order terms. 
Expanding in a power series in ?7,~2 V2 we obtain from 1)6.11) 

{e-\ = {e-''^),-n-HY,e-''^l + .... (6.2) 

Substitution of ()6.2j) in ()2.56j) yields for logp„ the expansion 

\ogpn = logpi°) + log(e"^0„ - (¥2)1 + . . . , (6.3) 

where for any function F of the angular variables we define the first order 
average 

(Fe-^^) 

The full average of a function F is given by 

(F) = {F\ - [ (FV2), - (F),(V2), ] + .... (6.5) 

If we can substantiate our claim that Vi is of order and n~^Y2 of higher 
order, then ()6.3|) and ()6.5|) are perturbation series in n~2. 

In this section we study V with the sole purpose of isolating an appropriate 
expression Vi. This will require an important amount of work. In section 
[7| we will show the anticipated result, viz. that for n 00 the average 
(^~^^)o tsnds to a finite constant C. In view of equation ()6.3|) this then 
establishes the asymptotic expansion of logp„ up to terms that vanish as 
n — >■ cxo. In sections [7| and |H1 we will determine first order expressions for 
various correlation functions. 
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6.2 An analog: the theory of elasticity 



Our determination of the interaction Vi will be based on the following con- 
sideration. The regular n-sided polygon is an obvious point of symmetry in 
phase space. In fact, Drouffe and Itzykson |27j constructed an exact lower 
bound for p„ by expanding the vertex vectors about their regular poly- 
gon positions. Such an expansion preserves, however, the long range order 
of the regular polygon and is therefore likely to cover only a small fraction 
of the effectively contributing phase space. Therefore we will proceed here 
as in the theory of elasticity, where one considers the deviations not of the 
positions, but of the interatomic distances, from their regular lattice values. 
The analogs of the interatomic distances are the angles C,m between the mid- 
point vectors and 7]^ between the vertex vectors. Hence we will expand V in 
powers of Sxm and 6ym defined in (14. 5|) . This procedure allows, in principle, 
for large deviations from the regular ra-sided polygon. In the nontrivial and 
rather lengthy analysis of this section we will step by step identify the order 
in n of the terms that we encounter, relegate to V2 those of order n~2 or 
smaller, and retain the others. The outcome of this will be expression ()6.39p 
for Vi. The terms that contribute to Vi are all quadratic in the 6xm and 

6.3 The interaction V as a series in the angles 

Definition ()2.54p of exp(— V) is a product of factors G', pm, T^^^^, and 
1 + n~^V. We will successively carry out the analysis of the orders in n 
appearing in each of these factors. Conceptually two different steps must be 
distinguished, even though in practice we have not strictly separated them. 
First, since the angles Pm,1m,C,m, and rjm are all small, even though not all 
of the same order, it will be convenient to expand all relevant expressions as 
power series in these angles. This is a matter of practical convenience rather 
than of principle. Second, we will convert the resulting series to an expansion 
in the 5xm and 6ym, retaining only the leading order. 

6.3.1 The factors G', Pm , and T^^~^ 

We first consider the function G. It is given by ()2.35|) in terms of the (3e and 
7£ which, in turn, are expressed by ()2.16|) in terms of Pi and of the and 
rjm- When differentiating G with respect to Pi at fixed ^ and r] we obtain 




m=l 




(6.6) 



m=l 
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where in the second step we used the geometrical sum rule ()2.12p to cancel 
the linear terms in the angles, and where the dots stand for terms of fifth 
and higher order in Pm and 7m- According to ()5.6p the sum + 7^ scales 
as ~ Since the sum on m contributes an extra factor n, we deduce from 
dnSD that 

G"(e,r/;/?,) = l + 0(n-i), (6.7) 

i.e., this derivative is to leading order independent of its arguments. The 
0{n~^) corrections in ()6.7|) may be treated perturbatively {i.e., they belong 
to V2), but will not be dealt with in the present work. 

Secondly, we obtain from ()2.37|1 and ()2.38|1 the expansion 

n n n 

n Pm = exp ( - J] + I J] + . . . ) . (6.8) 

m=l m=l m=l 

Expanding ()2.33|) and using ()2.37p yields for the series 

m m 

rm = r^ + \ Y^iPl - ll) - ^ - 7.') + • • • , (6.9) 

1=1 1=1 

which is valid for m = 1, 2, . . . , n; we recall that for m = n it is equivalent to 
the no-spiral constraint. 

Equation (j5.7|) tells us that (31 — 7^ and j3j — 7^ scale as ~ n~2 and n~2 , 
respectively. When m is of the order of n, as is typically the case, the sum 
on m contributes a factor n to the scaling, so that — t„ is to leading order 
proportional to n~2 . This leads us to define the auxiliary scaled variables 
by 

_i 

rm=n ^Tm, m = l,...,n, (6.10) 
which allow us to write the gauge condition as 

n 

5^r„ = 0. (6.11) 



m=l 



From (IT^ . (j237j), (ICTll . and (imj) we find that is related to the 
original midpoint distance Rm by 

Rm = -Rav(l + n 2rm) 

= i?, + 7r-^(rav + rj + 0(n-5). (6.12) 
In terms of scaled variables ()6.9|) becomes, if we use both equations ()5.7p . 

rm = rn + -^"^^ii^e - q) + y^a^AQ(^£ - q) + . . . , (6.13) 

e=i i=i 
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valid for m = 1, 2, . . . , n, and where the dots stand for higher order terms. 
Let us now consider ()6.8p . Since ~ and since the sum on n contributes 
an extra factor n, the first term in the exponent in ()6.8p is of order We 
will keep it. However, is of order n~2, so that its sum on m will vanish 
for n ^ oo and can be treated perturbatively. Hence the contribution to 
exp(— V) is a factor 



m=l 



m=l 



(6.14) 



Thirdly we consider the product of the T^^~^. Symmetrizing in j3m and 
7m with the aid of the no-spiral constraint G = we may write this product 

as 



sin(/3„ + 7 



\ iPm + 7m) COS f3m COS 7^ 
n 

E iUf^'m + ll) -W^ + lmf + ...). (6.15) 



m=l 



The sum Xlml-^m + 7m) contains n terms of order (n~^)^, is therefore itself 
of order n° and must be kept. Since /3m + 7m = im scales as n~^, the sum 
TlimiP-m + 7m)^ also scalcs as and may be treated perturbatively. The 
higher order terms in the sum in the exponential in ()6.15|) are at least quartic 
in the angles and may again be treated by perturbation theory. Hence, to 
leading order, 



m=l 



m=l 



(6.16) 



is the contribution of this product to exp(— V). 



6.3.2 The factor (1 +72- V)-" 

When writing ()2.39j) as we did, we anticipated that V thus defined would 
scale as ~ n° for n ^ oo. We will show here that this is indeed the case. We 
expand 

(1 + n-^Vy = exp ( - 1/ + in-V^ + . . . ) (6.17) 

and next wish to find \^ as a series in the angles. The decomposition ()2.39|) of 
A was motivated by the fact that in this way, and since the are quadratic 
in the angles, the series for V contains terms only from third order up. We 
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expand V as given by ()2.41|) and find 



m=l 

n n 

+ ~Y1 T-„(7„ + + — r^ilm + Pm+l)- (6.18) 

m=l m=l 

In the first line of (|6.18p the dominant term is the one where ^(7^ + Pm+i) 
is summed without any factor Tm- To estimate its order in n we may invoke 
fl5.5p : we will however perform a slightly refined analysis by using that 7^ + 
Pm+i = Vm, setting as before rjm = f] + Srjm with fj = 2Trn~^, and writing 

= lv{ll + Pl+i) + 'ISr]m{li + Pl+i)-bl^- (6-19) 

Because of the scaling ()5.H1 we have //(7m + Pm+i) ~ Hence, taking 

into account the effect of n times the sum on m in ()(i.l8|l . we conclude that 
the first term on the RHS of yields a contribution of order to V 

and should therefore be kept. The term — |?7j^ ()6.19|) contributes to order 

to V and can be treated perturbatively. The term |5?7m(7m + Pm+i) 
the RHS of ()6.19|) requires more attention. The factor Sr]^ is of the same 
order as fj in the first term, viz. ~ n^^. However, 61]^ is a zero average 
random variable and therefore |42j the sum on m will raise the order by only 
nz instead of n, whence it follows that the resulting contribution to V can 
be treated perturbatively. Putting these arguments together we get for the 
first line of I^H^ 

n n 

+ + + = + + ^(^"^)- (6-20) 

m=l m=l 

In the second line of equation ()6.18p we perform a similar analysis on the 
first sum. Using the sum rule (|2.38|) we can write 

n n 

- } , Tmilm + Pm+l) = - } Tmdrjm . (6.21) 
71 ^ — ' 71 ^ — ' 

m=l m=l 

Since Sr]m ~ ''7m ~ n'^, the product 5r]rn'Tm is of order n~^, and when summed 
on m with the prefactor w/tt in ()6.18|) . it would appear that it contributes 
a term of order na to V . This again overestimates the order. Since the 67]^ 
are zero average random variables, we expect jl7j in fact that the sum will 
raise the order by a factor of only , and hence that this term produces an 
order n° contribution to V. Therefore we keep this term. 
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In the second sum in the second hne of fl6.18|) we write T^(7m + Pm+i) = 
27m^^rl^ + T^^Srjm- Both these terms are of order (n~2)^ x = but the 
summation on m increases the order of the first one by n and of the second 
one, which is random jlTj, by only . Hence only the first term leads to a 
contribution of order to V and we get 

n n 

^ E + = E + ^(^"^)- (6.22) 

m=l m=l 

Finally, we must consider higher orders in V in the expansion ()6.17|) . 
Since we concluded that V is of order n°, it follows that all higher order 
terms n~^V^^, n~'^V^, ... in the expansion of the logarithm can be treated 
perturbatively. Hence the contributions to (1 + n~^V)~" that must be kept 
are given by (j6.2Up . (j6.2H) . and (j6.22j) . Upon adding them together we are 
led to 

n n n 

{l + n-'V)- = exp(-5^(7^ + /3^^,)-E^--^E^™^^™) 

m=l m=l m=l 

x(l + C(n"^)), (6.23) 
which completes our discussion of the factors composing exp(— V). 



6.3.3 Intermediate result for V 

Upon collecting the contributions identified above in equations ()6.7p . ()6.14|) . 
dEEl), and (jFT^ we find 

n n n 

m=l m=l m=l 

^ n 1 " 

= - E + + 2^-) + ^ E + ^(^"')' (6-24) 

m=l " m=l 

where in the second line we have passed to scaled variables. Although this 
expression has the appearance of a quadratic form, we recall that is itself 
a series in the angles starting out with a quadratic term. 

The important facts established at this point are (i) that V given by ()6.24p 
is of order n° as n — oo; and (ii) that the terms neglected in ()fj.24j) may be 
classified by their order in n~2 and can hence be arranged, if desired, in a 
perturbation series. 

A word of caution is in place concerning the averages (F). To arrive at 
V of equation ()6.24|) . we have used twice in section IB. 3. 21 that the sum of 
n sufficiently weakly correlated random variables has a random part which 
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is of relative order order n~2 with respect to its systematic part. For a 
two-point correlation at a distance k (measured in units of the midpoint 
index m) we expect that in an analogous calculation sums of k terms will 
intervene, for which the random part is of relative order /c~2 with respect 
to the systematic part. Hence the resulting expression ()6.24|) can be used 
for correlation functions at distances ~ n, but in the absence of further 
analysis not for k ^ . We speculate that in fact it is valid for correlations 
at all distances k rf for any e > 0, but our further work will not depend 
on that. 



6.4 Expressing V in terms of 5xjn and 5ym 

We return to the study of expression (j6.24|) for V. In order to prepare for 
integrating exp(— V) on the Xm and ym, as required by the average (exp(— V))^ 
in ()6.3p . we first wish to express the 6m,, Cm, and that appear in ()6.24|) in 
terms of the 5xm and 5ym- 

6.4.1 Recursion relation 

In principle, this should be done as follows. We may deduce from ()fj.l3j) the 
recursion relation 

rm-rm-i = rr^XmUm + 0{rr'^), m = l,...,n, (6.25) 

with ro = r„ as an additional condition expressing the no-spiral constraint 
and with the Um defined by 

Um = l{bm-Cm), 171 = 1,...,!%. (6.26) 

From ()6.26|) and ()5.4|) it follows that the Um satisfy 

Um - Um-i = n'^ fm-i , m = 2,...,n, (6.27) 

in which the fm are defined as 

fm-i = -l{5xm + 5xra-i-25ym-i), m = l,...,n (6.28) 

with the conventions Sxq = 6xn, Syo = and /o = /„. The extension of 
(I6.27|) to m = 1 defines uq. Because of the geometrical sum rules ()2.4|1 and 
(I2.1()|) we have mq = Un as well as 

n 

5^/m = 0. (6.29) 

m=l 

If initial values ro and uq are given, the system of equations ()6.25p and ()6.27p 
may be solved recursively in the midpoint index m, yielding the and Um 
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in terms of the fm, and hence in terms of the Xm and i/m- In the end the 
initial values are determined by two conditions, viz. the no-spiral constraint 
r„ = ro and the gauge condition ^^^^ = 0. Proceeding this way gives a 
host of terms many of which contain multiple sums on the midpoint index. 
Their order in n may be determined by the methods of section 16.31 This 
route, however, is cumbersome as too much detailed information is carried 
along. We will therefore adopt a different but equivalent approach which 
selects the leading order terms more easily. It also reestablishes the cyclic 
invariance in the index m and the symmetry under reversal of orientation. 

In fl6.25|l we substitute Xm = x + Sxm- Of course 6xm is of the same order 
as X. But since the 6xm are independent zero average random variables, by 
a reasoning analogous to the one advanced in subsection ()6.3.2|) . the random 
part, when summed on m [as for exemple in ()6.13|l ]. is of order relative 
to the systematic part. It may therefore be neglected to lowest order. Upon 
combining ()6.25|) and ()6.27|) and setting x = 27r we then find for the 
second order recursion 

3 

Tm+l - 2rm + Tm-l = '2Tm~2 f^^ (6.30) 
which holds the key to much of what follows. 



6.4.2 Solution in terms of Fourier transforms 

In order to solve the recursion (16.301) we define the Fourier transforms 



r^ = n-iJ2 e^'^'^^^/^r^ , Ug = n-^ ^ e^-^''^'"-^)/"^^ , (6.31) 



m=l m=l 

and 

1 " 

/, = ^ V e^-'^"^/"/^ , (6.32) 

ra=l 

where g = 0, ±1, ±2, . . . , ±(|n — |) for n odd and g = 0, ±1, ±2, . . . , ±(|n — 
l),|n for n even. Here q is the 'wavenumber' and 27r/g the 'wavelength' 
along the cell perimeter, the scaled index 27rm/n being equal to the polar 
angle in the plane up to corrections of order In Fourier language the 

unknowns Um and of (j6.27p and (j6.30|) are readily solved in terms of the 
forces fm and one finds 



Uq = --. ^^fq - --fq, ^7^0, (6.33) 

m sm — iq 



vr^ , 1 



r. 



^ --of,. 5^0. (6.34) 



sm ^ q'^ 
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The no-spiral constraint and the gauge condition become uq = and fo = 0, 
respectively. The ~ signs in ()6.33p and ()6.34|) indicate the limit n ^ oo 
at fixed q and amount to neglecting higher orders in q/n; these are small if 
we restrict ourselves to q values on a scale ~ n^"*^, i.e. to spatial distances 
which, in units of the midpoint index m, scale as ~ n'^ with e > 0. 
Next, Fourier transforms of Sxg and Syg will be defined as 

^ n n 

Xg = T "E e'"''""/"^^^ , Yg = V e^-Wn^^^ . (6.35) 

From and (lOKl) it then follows that 

/, = -i(l + e-2-'''/")X, + F, ^ -(X,-y,). (6.36) 

We remark that the factor n~5 in ()6.35|) is the usual one for a Fourier 
transform of n independent variables; it guarantees that Xg and Yg are of 
order n°. The factor in the Fourier transforms ()6.31|) and ()6.32|) then 
follows from the requirement that in each of the equations ()6.33j) . ()6.34|1 . and 
fl6.32|l all members be of the same order in n. 

6.4.3 Defining Vi 

Employing ()6.3H) and ()6.35|) we are able to express V of equation ()6.24|) as 

Y = Y^ M-. + 2r,r_, + 2f,yiJ + 0(n-^), (6.37) 

where q runs through all nonzero integers and where we Fourier transformed 
using that and — are both equal to Um + 0{n~2). We now insert the 
solutions ()6.33p and (j6.34|) in ()6.37|) . We relegate to the term n^^y^ in ()6.ip 
also the higher orders in q/n discussed above and arrive at 

V = Vi + 0(r2-^) (6.38) 

where we finally decide to define Vi as 

In fl6.39|) the sum on q converges due to the presence of the factors 1/q^ and 
This confirms again that Vi is of order n°. The contributions to Vi 
come essentially from wavelengths at the scale of the cell itself and are fully 
separated in scale from the independent local degrees of freedom that entered 
the noninteracting calculation. 
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In this section we have extracted from V, defined by fl2.54j) . its leading 
order term Vi, given by ()6.39|) . With this expression at our disposal we can 
now go ahead and calculate the consequences of the interaction. Expression 
(I6.39|) could be considered as an elastic Hamiltonian in the harmonic ap- 
proximation, if it were not for the fact that the Xg and Yg are non-Gaussian 
variables. It will appear below that the deviations from Gaussian statistics 
affect the results of interest only to higher order in n~ 2 . 

7 First order 

At the end of section IHl we concluded that V is of order nP. This means, 
more precisely, that it is ~ in the region of phase space anticipated to 
contribute dominantly to the average (exp(— Vi))^. We therefore expect this 
average also to be of order By showing that below we will accomplish 
an essential step. The development of this and the following section are 
standard mathematical physics. 

7.1 Calculation of pn 

We consider the quadratic form Vi given by ()6.39|) above. It is useful to set 

Zg = iXg,Yg) (7.1) 

and to write ()(i.39|l more compactly as 

Vi = 5^Z,-V,-Z!^, (7.2) 

where the superscript T indicates transposition and where 

The exponential exp(— Vi) with Vi given by ()7.2|) has to be substituted for 
exp[—W{x,y)] in ()3.5|) . The variables of integration Xm and ym are now 
coupled by Vi. In order to decouple them we apply a Hubbard-Stratonovich 
transformation. We write 

V',= (i^i^) (7.4) 
and represent exp(— Vi) as the integral 

(detV. 



q^O q^O q^O 

(7.5) 
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We have that det Vg = —q ^ is negative and the paths of integration of the 
if)^ and tp^ in the complex plane should be chosen such that the integrals in 
()7.5|) exists. Let and ip^ be related to tp^ and t/'^ in the same way that 



5xm and are related to Xq and Fq, i.e. by ()6.35|) . We set ipQ = ipQ 



0. 



Then 



J2 {r^x^q + i^iY^q) 



(2vr) 



1 " 



(7.6) 



m=l 



We insert fl7.6|) in ()7.5p and the integral representation ()7.5|) in ()3.5|) . We 
commute the integrations to the outside and are then able to perform the 
integrals on the Xm and i/m- The result is 



Pv 



(detV-i)i 



X 



2n 



exp 



n 

L27ri 



„2nnu " 



dw 



2n 



n 1 



2; \ -2 
m 



m=l 



n 



zvri 



n 1 



m=l 



(7.7) 



Hence one sees that the u and f integrals factorize again. The products on 
m represent an insertion of the type h{u, v) discussed in section 14. H which 
may be taken at the saddle points u = Uc = tt"^ and v = Vc = {2tx)~^. Doing 
so and exponentiating these products we get 



Pr. 



(detV^ 



An 



X 



n 

exp[-l5^^,>^^-5^27^(C 

gj^O m=l 



X exp 



g^O 



(7.8) 



m=l 



where pn^ was given in ()4.1|) and the 'remainder' TZ{x) stands for the infinite 
series 



^ 1 • I 



fc=3 



(7.9) 



We will delay to section 17.2.21 the proof that in the limit n ^ 00 the remain- 
ders in (17. 8|) contribute only to relative order n~^. Admitting this here and 
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setting 



we therefore have 



E 



An 



(detV-^)^ 
^^0 [det(E + V^-i 



1 
2 



11 dV', 



exp 



(7.10) 



iJ]V',-(E + v-i)-V'. 



9=1 

where we abbreviated 



(7.11) 



det(l + V,E) 



(7.12) 



with 1 standing for the 2x2 unit matrix. Upon combining (|7.11|) and (|7.12j) 
with expression ()4.ip for one finds for Pn the result announced in ()1.2|) - 
(jl.3j) of the introduction, 



C (8 



2\n 



47r2 (2n)! 



■[1 + 0{n-'^)]., 



where 



9=1 



(7.13) 



(7.14) 



In the course of the calculations of the previous sections we have encountered 
correction terms of relative order n~2. In equation ()7.13j) we have included 
these 0(n~2) terms again, although strictly speaking we have not shown that 
their net effect does not happen to vanish. 

This is the first analytic determination of the large n behavior of p„. 
Conjectures in earlier work (see [^) have almost invariably had the form of 
decaying exponentials, whether or not multiplied by a power of n, or sums of 
such expressions. Only Drouffe and Itzykson |2ZI come close to our analytic 
expression: they assume that to leading order logp„ ~ an and estimate 
a ^ — 2, a value confirmed by (j7.13p . We refer to for a plot of the 
asymptotic expression ()7.13|) and a comparison with Drouffe and Itzykson's 
simulation data. 
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7.2 Correlations to first order 



The calculational technique that led to the expression for p„ in the preceding 
subsection, may also be used to find correlation functions at long distances. 

7.2.1 Long distance correlations along the perimeter 

We will study the correlations between the angles and rjrn in terms of their 
scaled counterparts xi and i/m- Passing to Fourier transforms with the aid of 
f|6.35|) we have 

(SxeSx^) = ^ e2-^(^-^)^/"(X,X_,) (7.15) 

and analogous relations for {SyiSym) and the cross correlation {Sxi6ym) ■ From 
definition ()6.4|) and equations ()7.5p and ()7.6|) it is clear that we obtain the 
first order expression by inserting in the integrand of (j7.5|) a sec- 

ond derivative —d'^/dtp^dip^^ that acts on the factor exp(— i -^^-Z^g). 
Carrying the calculation through we find that we have to differentiate the 
exponential factors in the second line of ()7.8|) . whence we have, for g 7^ 0, 

Here {■■■)(. is a Gaussian average with respect to the weight est x exp [ — 

\Yliq^o'^q' ^ + V^^)-'0_J and the dots represent a well-defined series of 
quartic and higher Gaussian averages coming from the expansion of the re- 
mainder TZ in (j7.8p . When evaluated, these higher order terms contain factors 
n~^, as will be shown in section 17.2.21 After completing the calculation for 
the other basic correlation functions we obtain 

{yqy-q\ ^ 1 - {rqr-q)a , 

{X,Y^q\ ^ - -Mr--q)a^ (17^0^ (7.17) 

where the ~ sign indicates, as before, asymptotic equality for n —>■ 00. In 
view of ()fj.5j] we may drop the index 1 on the averages in ()7.17j) . The three 
Gaussian averages in ()7.17|1 are due to the interaction Vi. They come in 
addition to the constants (|, 1, and 0) that represent the Fourier transformed 
noninteracting averages ()4.6|) . Straightforward calculation yields 

{rqr-q)a = -4A-'g-^ ^t^o, (7. is) 
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From ()6.5|) . ()7.18p . and ()7.17p we obtain for the basic correlation functions 
the following set of first order expressions, 

^ l + 2A-i(g-2-g-4), 

^ 2A-lg-^ g^O, (7.19) 

which represent the main result of this subsection. The terms with the factor 
in ()7.19p have their origin in the interaction Vi. The important point 

about the ipip correlations is that they all decay at least as ~ when 
g — s> oo. This will turn out to be an essential element in our analysis below. 
Equations ()7.19|) may be inverse Fourier transformed according to ()7.15|1 and 
its analogs so as to yield expressions for the spatial correlations {6xe6xm), 
{6ye6ym), and {5x^5ym)- We will work these out in detail only for the more 
restricted correlation function {{Xg — Fg)(X_g — Y^q)) = {fqf-q) to be con- 
sidered in section |H1 Because of our discussions on spatial scales at the end 
of section IB. 3. 31 and following ()6.33|) - ()6.34p . the real space correlations so ob- 
tained are valid only at distances large compared to the individual perimeter 
segments. 



7.2.2 The remainders TZ 

Having gained control of the correlation functions, we can now prove our 
claim that the remainders TZ in the exponentials in ()7.8p contribute only to 
higher order in n~^. We may expand these exponentials as perturbation se- 
ries in the ip^ and ip^^. The odd terms vanish because of symmetry. The 
lowest order contributing terms are therefore proportional to '^^i^'m)^ 
J2mi'^rn)'^- Let US cousider the latter. It leads to a correction term which, 
relative to the term calculated, is of order ^m((V'm)'^)G- therefore eval- 
uate 



n 

where we used Wick's theorem. The fact that {''pq^'-q)^ decays as q^^ is cru- 
cial here. It means that the sum on q inside the brackets in ()7.20j) converges, 
and hence ^m((V'm)'^)G order n~^. We therefore had the right to neglect 
the TZ terms in ()7.8|1 . The other terms in the perturbation series in the 
and can be argued in a similar way to be at least of order n~^. 
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7.2.3 Debye-Hiickel 



We remark that in order to arrive at the above results, as well as those that 
will follow in section |H1 one may bypass the Hubbard-Stratonovich transfor- 
mation by considering the variables of integration Xq and Yq as zero mean 
Gaussian variables of variances determined by those of {5x1^^ and {Sy'^)^ 
given in ()4.2|) . i.e. with weight 

exp {-J2^q-E-'-Z^q) = exp (- 5^X,X„, - (7.21) 

The procedure of replacing non-Gaussian variables by Gaussian ones of the 
same mean and variance is known in plasma physics as the Debye-Hiickel 
approximation; it is exact in the limit of high temperature. This limit hap- 
pens to apply to our case, where the role of the temperature is played by 
n: the interaction 'energy' Vi would be extensive {i.e. proportional to n) 
in an ordinary thermodynamical system, but is here of order n°, which is 
equivalent to having an inverse-temperature prefactor equal to n~^. 

In the preceding subsection we have chosen to go through the full Hubbard- 
Stratonovich transformation. This has allowed us to show that the procedure 
is exact and to prepare the way for the calculation of higher order terms, if 
desired, in the power series expansion in n. 

8 Weight functional for the perimeter 
8.1 Weight function for the fq 

Given that for large n the perimeter approaches a circle, it is natural to ask 
what the probability distribution is for its deviations from circularity. This 
question has prompted work, in particular, in the mathematical literature 
inni EZj ; but within the framework of large-cell limits different from ours. 
Deviations from circularity may be characterized in various ways, e.g. as 
the difference between the radii of the smallest circumscribed and the largest 
inscribed circle [5^ IH7j . Here we will exhibit in full generality the weight 
of a specific trajectory of the perimeter, characterized by the set of scaled 
'excess midpoint distances' {ri, r2, . . . , r„}. The Fourier transforms fq of 
these distances are related to the differences Xq — Yq by q^fq = — fq = Xq — 
Yq., which follows from (j6.34j) and (|6.36|) . Hence one of the variables Xq and 
Yq is redundant for the description at which we are aiming. The probability 
distribution of the Xq and Yq is, up to normalization, equal to the weight 
exp(— Vi) with Vi given by (j6.39p . In order to derive from it the marginal 
distribution Pi{{fq}) for the /g, we short-circuit the Hubbard-Stratonovitch 
procedure in the way indicated in section 17.2.31 Upon integrating exp(— Vi) 
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directly with the appropriate noninteracting Gaussian weight ()7.2H) we find 
the first order probabihty law 

Pi({/J) ^ est X exp (^-lJ2^JJ-g) (8.1) 

with Aq given by (j7.12j) and where, because of ()(i.29|l . it is understood that 
/o = 0. The law (jH.lll leads to the correlation 

(/J,')-|A-\+,',o, q,q'^0, (8.2) 

which may alternatively be derived from ()7.19j) . We remark parenthetically 
that the weight 1)8.11) is suitable for computer simulation of the large scale 
structure of Voronoi cells. The rest of this section elaborates on equations 
(lO) and dHU. 

8.2 Differential equation for R{(j)) 
For 72 oo it is possible on those scales to consider the angle 

27rm 



n 

continuous variable ESI and to write 



(8.3) 



rm = ri^), fm = ^fi^). (8.4) 

Recursion ()6.30p relating the to the fm then becomes the second order 
differential equation 

r"(0) = /(0). (8.5) 

Here /(0) is Gaussian noise whose correlation function (/(0)/(0')) "^i^^ 
now determine. To this end we set 0' = 2nm'n~^ and consider {fmfm'), 
which with the aid of ()6.32|) may be written as 

{fmfm') = ^J]J]e-^-'('^-+'''-'V"(/J,,) 

n 

q q' 

- ^E^"'^'"''^"' (8-6) 

where to arrive at the second line we used ()8.2p . Passing from to /(0) 
and using (|7.12p we obtain from (|8.6p 



(/(0)/(0')) ^ |[5^e-'''(^-^') + r(0-0') 

q^O 

= |[27r5(0-0')-l + r(0-0')], (8.7) 
44 



where the 27r-periodic function r(0) is given by 



oo —2/1 —4 

r(0) = 2 5^^^-;^cos#. (8.8) 

9=1 

This correlation function resuh merits a comment. It is instructive to cal- 
culate, for comparison to ()8.6p and ()8.7p . also the simpler noninteracting 
average (/m/m')o- To this end one expresses fm and fm' in terms of the 6x£ 
and 6ye using ()(i.28j) and employs equations ()4.fi|l . In the continuum limit it 
then appears that in ()8.7|1 the terms 2TT6{(j) — 0') — 1 represent the nonin- 
teracting result and r(0 — </)') the contribution generated by the interaction. 
The expression for this contribution was derived under the restrictive con- 
dition that 10 — 0'| ^ 27in~^. We speculate that ()8.7p is vahd, in fact, for 
10 — 0'| scaling as ~ n~^^'^ with arbitrary e > 0. Finally we remark that 
expressions similar to (j8.7p for the correlations ((5x(0)(5x(0')), {6x{(f))6y{(j)')), 
and {6y{(f))6y{(l)')) may be found in a fully analogous way. 

In the continuum limit the Fourier transform of Ag becomes 

J2 e^'^^^-^'^Ag = 27r(5(0 - 0') - 1 + f/(0 - 0') (8.9) 
with the 27r-periodic function t/(0) given by 

oo 

f/(0) = -2^(g-2 -4g-^)cosg0. (8.10) 

q=l 

Similarly, in (j8.9p the term f/(0 — 0') represents the effect of the interaction. 
Both r(0) and U{(f)) are of zero average on the interval [0, 27r]. 

The sums on q in equations ()8.8|) and ()8.10|) may be carried out analyt- 
ically [see reference jSU], equations (1.443.3), (1.443.6), and (1.445.2)]. One 
finds that ()8.10p leads to the quartic polynomial 

t/(0) = -|7r2(^7^2-l) + |(K-l)(0-7^)^-|(0-7r)^ < < 2vr, (8.11) 

and its 27r-periodic continuation. We have written it such as to bring out its 
symmetry under the replacement i— > 27r — 0. The potential U{(f)) in ()8.10p 
consists of two terms which merit separate discussion. The the first term is 
the one-dimensional electrostatic interaction but with opposite sign, so that 
it is attractive for like charges; and the second one is repulsive. Near the 
origin we have the expansion 

f/(0) = f/(0) + 7r|0| - (1 + |7r2)02 + . . . , (8.12) 

in which the linear term comes from the electrostatic interaction. As a result 
[/(0) has maxima at = ±0^, where 0^, = vrfl — {1 — 3/(27r^)}2] = 0.079...7r. 
We will not exhibit the analytic expression of r(0) but have represented its 
graph in figure El It has minima at = ±0j,, where 0j, ~ 0.139...7r. 
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Figure 3: The functions and T{4>) of equations (|8.11|) and H8.8|) . respectively. 

8.3 Weight functional for 

Using ()8.9|) we subject the quadratic form in the exponential of ()8.H) to 
an inverse Fourier transformation. In the continuum limit the probability 
distribution -Pi({/q}) becomes a functional of /(0), or, because of ()8.5|) . of 
r"(0). We find 

Pi[r(0)] ~ est X exp ( - i/[r"(0)]), (8.13) 

DVT Jo 

+ A^i r#2r"((/)i)f/(</>i-</>2)r"(</.2). (8.14) 

l^TT^ Jo Jo 

Here the last line, which is the contribution from the interaction V, is anal- 
ogous to the Coulomb energy of a line charge density r"(0) in a periodic 
one-dimensional system. Equations ()8.13p and ()8.14|) are valid for functions 
r(0) coarse grained to a scale much larger than the distance between the 
individual vertices, i.e., to angle differences ^ 27rn~^. Moreover, the 
sum rules encountered during its derivation now lead to three restrictions on 
the functions r(0) that enter ()8.14|) : 

(i) The geometrical sum rule ()6.29|) . taken in the continuum limit and 
combined with the differential equation ()8.5p . implies the 'charge neutrality' 
condition J^^ dcj) r" ((p) = r'(27r) — r'(0) = 0. It follows that we may add an 
arbitrary constant to the potential U in (j8.14p . 
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(ii) The no-spiral constraint becomes /o d0r'(0) = r{2n) - r(0) = 0. 

(iii) The gauge condition takes the form Jg'^d0r(0) = 0. 

This completes our discussion of the first order probability functional Pi. 

We now wish to reconvert the above findings to results for the original 
midpoint distances Rm = R{<P) = Rc + SR{(f)). From ()6.12|) taken in the 
continuum limit and combined with (jl.lUp we find 

6R{(f)) = 7r-5 [r,, + r(0)] + (9(^-5). (8.15) 

It follows that the probability distribution of R{(f)) is determined by the 
product of Pav('^av) givcu in (j2.45p and Pi[r(0)] given in (|8.13p - (|8.14p . Still 
setting 

F(0) = 7r"^/(0) (8.16) 

we are led to the results ()l.lH) - p.l3|) announced in the introduction, in which 
the factors 4A have been restored. 

8.4 An application: fluctuations of Rm around i?av 

As an application we determine how, within a single realization of an n-sided 
cell, the midpoint distances Rm deviate from their average i?av This extends 
the considerations of section B. 4. 41 where we showed that R^v is itself within 
order n° from its ensemble average Rc = ir'^n^ . 
We have straightforwardly 

n n 

n-' Y.{{Rm ~ R..r) = n-'RlY.^rl) 

m=l m=l 

= n''RlY.^~'{Ul,) 

q OO 

= + (8.17) 

9=1 

where the first step is based on the first line of ()6.12|) . in the third step we 
used ()(i.34|l . and in the last step we inserted the explicit expressions for Pav 
and {fqf^q). Equation ()8.17|1 yields the result announced in p.9j) . 

8.5 Convergence of R{(f)) and S{(f)) 

The set of midpoints {Rm \ m = 1, . . . , n} and the set of vertices {S^ | m = 
1, . . . , n} each define an n-sided polygon. In the limit of large n these poly- 
gons become the curves -R(0) and 5(0), respectively. Of these, S^cj)) is the 
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true perimeter. We show here that for n — oo the curve R{(f)) converges 
rapidly to S{(j)) so that in the discussion above we were justified in referring 
to R{(f)) as the perimeter. 

Two points Rm = {Rm, ^m) and = (Sm, ^m) have angular coordinates 
such that |$m — = 1 7m I = 0{n^2^, Since increases by steps of order 
~ when i increases by 1, there exists a third point R^' = {Rm',^m'), 
with m' depending on m, such that 

= C(n-^). (8.18) 

Typically we will have \m' — m| = 0{n^). In order to estimate how the 
distance between R^' and Sm scales with n we write 

|R„, - S„|2 ~ - S^Y + Rli<^^m' - ^mf, (8.19) 

where we used that to leading order Rm' — Sm — Rc- In order to deal with 
the first term on the RHS of ()8.19|) we consider the triangle inequality 

Sm\ ^ \Rm' Rm\ ~l" \Rm Sm\- (8.20) 

Since Rm = Sm cos 7m, we have 

\Rm - Sm\ ^ \Rcll = 0{n-^2). (8.21) 
Furthermore, from ()6.25p it follows that 

m' 

rm' - rm - 27m~^ ^ m 

e=m+i 

~ 27m"^\m' — m\um , (8.22) 

where in the last step we used that for i in an interval of length \m' — m\ ~ na 
the value of U£ varies negligibly with respect to its typical value. Hence ()8.22p 
leads to \rm' — rm\ = 0{n~2)^ which together with (j6.12|) gives 

\Rm' -Rm\ = Rln-"^ \rm' " r^l = 0(n-5). (8.23) 

Upon using and (j^T^ in (j^T^ we find that \Rm' - Sm\ < 0{n-^). 

When substituting this in (j8.19p and using (j8.18p one obtains that |Rm' — 
Sm| = 0{n~^). This says that any vertex on the polygon through the Sm 
is, typically, within a distance n~2 of a vertex of the polygon through the 
Rm- In the limit — > oo we may therefore indifferently refer to either of 
them as the perimeter. 
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9 Lewis' law 



The average area A of a. Voronoi cell is equal to A = A~^, where, here and 
henceforth, we restore the particle density A. Its average number of sides (n) 
is topologically constrained to (n) = 6. One may now investigate the average 
area An of an n-sided Voronoi cell. Lewis' empirical law 2\ states that An is 
linear in n. The law is often expressed with a single free parameter c. 

An = [c{n ~ 6) + l]X-\ (9.1) 

so as to imply that Aq = A^^, i.e., the average area of a cell having (n) 
sides is equal to the average area of all cells; this appealing identity is not, 
however, a mathematical necessity. The law ()9.1|) accounts very well for data 
on cell areas coming from a wide range of different experimental situations 
|3j, despite deviations often found for the small-n values n = 3 and n = 4. 
Lewis' law also seems to apply j2Sl HH IHH lEB] to models different from 
Poisson- Voronoi diagrams, whether dynamical or not. 

The numerical estimates of the coefficient c presented in reference jH] are 
in the range 0.199 — 0.257. In particular, Lauritsen et al. |28j estimate c = 
0.23; Drouffe and Itzykson [22], on the basis of their Monte Carlo simulations 
for large n, observe that c is 'very close' to |; and Miles and Maillardet (Mj . 
on the basis of a comparison of the cell area with the area of its fundamental 
domain |40 , conjecture that c 'does not differ greatly' from |. It has been 
noticed [38^ that when the 'natural' boundary condition ^2 = is imposed 
on ()9.H) . one finds by coincidence c = | exactly, whence 

An = \{n-2)\-\ (9.2) 

Various 'derivations' of Lewis' law (j9.1|) have been proposed jSHl IHH 1^ - 
mostly involving the maximization of an entropy defined as a functional 
of the probability distribution {pj | j = 3, 4, . . .}. None of these derivations, 
however, is based on first principles, and Chiu [HHIIS^ shown that entropy 
arguments of this kind do not allow firm conclusions. Moreover, numerical 
data on ^3 and ^4 for Poisson- Voronoi diagrams clearly show deviations 
from linearity and hence rule out laws such as ()9.1|) and ()9.2j) for the full 
(nonasymptotic) regime. 

It suffices to extend the arguments of the preceding sections only very 
slightly to arrive at a new analytic result. Let us look at the midpoint 
distances Rm- By combining equations (j2.29j) . (j2.37j) . and (|6.10|) we find 
that Rm = -Rav[l + C(?^~2)]. Furthermore the discussion of subsection 12.4.31 
has shown that R^^ = Rc[l + 0{n-^)]. It follows that Rm = Rc[^ + 0{n'^)] 
for all m, where we recall that Rc = [n/A-KX]^. Hence the cell perimeter 
stays within a distance of order unity from a circle of radius ~ na. Therefore 
the area of the n-sided Voronoi cell is sharply peaked around an average An 
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for which we find the asymptotic relation 

^ TiRl = -. (9.3) 

Hence we have demonstrated from first principles that Lewis' law holds to 
leading order in an asymptotic expansion for large n, and that its coefficient 
is c = i. 

10 Comments and conclusion 

By setting out to find the probability p„ for a Voronoi cell to have n sides, 
we ended up unraveling the detailed statistical structure of the many-sided 
Voronoi cell. We have obtained the asymptotic expansion of logp^ up to 
terms that vanish as ri — > oo, and have shown that a perturbation series 
may be set up for the higher order terms. We found that the problem is 
characterized by two spatial scales which in the limit n — oo, and to the 
order that we calculated, become completely disjoint. The first one is the 
microscopic scale and the second one the scale of the system size. 

The exact methods and results of this work suggest, in retrospect, several 
heuristic considerations. These may be expected to lend to our arguments 
a fiexibility that will allow them to be applied in more general situations. 
We reserve the heuristics, together with further extensions of the present 
results, to a companion paper [32]. Our work also opens the way to efficient 
Monte Carlo simulations of large Voronoi cells; we also leave this and other 
numerical applications for future work. 

Voronoi cells in spatial dimensions higher than two have likewise received 
a great deal of attention in the literature jSj . We do not see how to calculate 
an analytic result analogous to our equations ()1.2j) - ()1.3|) in dimensions above 
two. The reason is that the idea of constructing a Markov process along the 
cell perimeter, which was at the root of the present work, has no obvious 
extension to higher dimensions. Nevertheless, we believe that it is possible 
without detailed calculation to transpose many of the scaling relations found 
in this work to higher dimensions. 

The problem of random points in a plane considered here (the Poisson 
point process) is twin to that of random lines in a plane: the Poisson line 
process. Instead of Voronoi cells, random lines define Crofton cells, about 
whose statistics which one may ask similar questions. In recent mathematical 
literature Voronoi and Crofton cells have often been discussed together 

EZ] • In physics random lines in a plane were was introduced by Goudsmit 
|58j . who was interested in the statistics of tracks in a cloud-chamber. He 
asked, specifically, about the probability that three or more independent 
tracks intersect at almost the same point and converted this question into 
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one about the area of the Crofton cell. We are confident that the methods 
of this paper can be made to work for the description of the large n-sided 
Crofton cell. 

In metallography large cells were observed by Aboav 29^ in the arrange- 
ment of grains in a polycrystal. Aboav's empirical law (also called the Aboav- 
Weaire law) states that the neighbor of an n-sided cell has itself on average 
'^n = flo + o,in~^ neighbors, where oq and ai are constants. This law has 
been widely applied to most of the experimental and theoretical systems dis- 
cussed in the introduction. The law is nevertheless violated, although only 
slightly, by numerical results on Voronoi diagrams On the basis of 

the analysis of the present work we will in a future study [30] find the true 
expression for m„ when n is asymptotically large. 
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